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ABSTRACT 

We explore the gravitational influence of pressure supported stellar systems on the 
internal density distribution of a gaseous environment. We conclude that compact 
massive star clusters with masses > 10 6 Mq act as cloud condensation nuclei and 
are able to accrete gas recurrently from a warm interstellar medium which may cause 
further star formation events and account for multiple stellar populations in the most 
massive globular and nuclear star clusters. The same analytical arguments can be used 
to decide whether an arbitrary spherical stellar system is able to keep warm or hot 
interstellar material or not. These mass thresholds coincide with transition masses 
between pressure supported galaxies of different morphological types. 

Key words: galaxies: star clusters — galaxies: fundamental parameters — globu- 
lar clusters: general — globular clusters: individual: uj Cen — ISM: general — ISM: 
kinematics and dynamics 



1 INTRODUCTION 

Stars are believed to form during one single event in com- 
pact star clusters with a total mass ranging from a few solar 
masses to a few million solar masses. But in recent years ob- 
servations have revealed that stellar populations of compact 
star clusters are more complex than a single aged stellar 
ensemble. 

V ery young star clusters such as th e ONC (|Palla et al.l 
2005) and a Orionis (|Sacco et al.ll2007l ) may probably har- 
bour a few low mass stars which are more than 10 Myr older 
than the main bulk of the ir stars. This may be the result of 
extended star formation (|Palla et al.l 120051 ) or the capture 
of stars, which are born in surrounding former star forma- 
tion eve nts, through the deepening potentia l during cloud 
collapse ( Pfla mm- Altenburg fc Kroupall2007l ). 

On the other hand some of the old massive globular 
clusters of both the Milky Way and Andromeda exhibit a 
spread in metallicity and/or have subpopulations with a he- 
lium overabundance suggesting the occurrence of multiple 
star formation events in these compact star clusters in the 
past (see Section [2] for a compilation of observational evi- 
dences). 

All star clusters with evidences for multiple stellar pop- 
ulations have in common that their total mass is around 
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10 Mq or higher. Furthermore, star clusters with a mass 
larger than ~10 6 Mq show a correlation between the cluster 
mass and the cluster size, whereas star clusters less massive 
than ~10 6 Mq have a constant radius dWalcher et al]|2005l : 
lHasegan etal]|2005l : Ibabringhausen et al.ll2008l ). 

Mass thresholds of globular clusters have been 



considered by several authors (Morgan fc Lake! 1 19891 : 
IShustov fc Wiebd l2000l : iRecchi fc Danzigerl |2005| )~ address- 
ing the question how massive a star cluster must be in order 
to keep ejected material by AGB stars or supernovae. How- 
ever, these thresholds can not explain the su ggested growth 
of nu clear star clusters by repeated accretion (jWalcher et al.l 
and large internal spreads in metallicity. 

iLin fc Murravl |2007l ) discussed the collective accretion 
of gas by globular clusters onto their member stars embed- 
ded in a 10 6 K hot back ground gas corresponding to todays 
circumstances. But globular clusters are assumed to have 
formed in denser environments, as for example nuclear star 
clusters are fo r ming today. Additionally, the treatment by 
ILin fc Murravl (|2007l ) does not reveal the existence of a mass 
threshold at 10 6 Mq. 

In this paper we show that compact massive star clus- 
ters with masses > 10 6 Mq are able to accrete gas from an 
embedding warm inter stellar medium. We first summarise 
the observational evidences for multiple stellar populations 
in massive star clusters in Section[2] The theoretical criterion 
for possible gas accretion by massive star clusters is derived 
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in Section [3] and the 10 6 M© accretion mass threshold is 
obtained in Section [4] Additional mass thresholds are found 
and identified with transitions in the radius-mass plane of 
spherical and pressure supported stellar systems (Section[B]). 



2 EVIDENCE FOR MULTIPLE POPULATIONS 

The most massi ve globular cluster, ui Cen, with a mass 



of ~2.5-10 6 M Q ijMcLaughlin fc van der Marelll2005l) show s 



a wide spread in meta llicity jFreeman fc Rodger* 1 19751) . 
iHilker fc Richtlerl |2000) determined by Stromgren photom- 
etry three different stellar populations being ^1-3 Gyr and 
~6 Gyr younger than the oldest one. iBedin et alj (|2004T ) 
confirmed the existence of multiple evolutionary sequences 
in W Cen. 

The properties of RR Lyrae stars observed in 
the globular clusters NGC 6388 and NGC 6441 with 
masses of ~l.lxl0 6 Mr and ~1.6 xl0 6 M©, respectively 
|McLaughlin fc van der Marelll2005h . can b e reproduced by 
the composition o f two distinct populations l|Ree et al.ll2002i ; 
lYoon et alJliooi). 

Recently, iMilone et all (|2007l ) showed that the split 
of the sub-giant branch of the Galactic globular 
cluster NGC 1851 with a ma ss of ~3.1xl0 5 M© 
ijMcLaughlin fc van der Marell l2005h corresponds to the ex- 
istence of two different stellar popu lations with a n age 
spread of about 1 Gyr. Furthermore, IMilone et all (|2007l ) 
concluded that the observed RR Lyrea gap in NGC 1851 
req uires an age dif fe rence of ~2-3 Gyr. 

iLehnert et all l|l99ll ) reported a possible metal- 
licity spread in M22 with a mass of 4.4- 10 5 M, 



McLaughlin fc van der Marej|2005h , whereas iRichter et al.l 



1999T I found no evidence for a metallicity spread using 
Stromgren photometry. 

Metallicity spr eads are also reported for Gl 
l|Mevlan et al. l200l[ ). the most massi ve globular clus- 
ter in M31 wit h a mass of ^8T0 6 Mm |Baumgardt et al.l 
120031 ') and M5 4 JSaraiedini fc Lavdenlll995l) with a mass of 
~2.0-10 6 M© |McLaughlin fc van der Marell [2005h . 

Evidence for metallicity spreads in the three massive 
glo bular clusters G78, G2 1 3 and G280 of M31 are reported 
by iFuentes-Carrera et al.l (|2008l ). These clusters have in- 
terna l velocity dispersions as high as Gl l|Diorgovski et all 
1 19971 ) and therefore must have masses larger than 10 6 M©. 

Various explanations for the origin of multiple stellar 
populations in globular clusters and abundance anomalies 
such as ejecta from AGB stars (eg. IVentura fc DAntona 
2008) or r otating mass ive stars (eg. iDecressin et al.l 120071 ) 
exist fsee lRenzinill200sl for a summary). Self-enrichment by 
rotating massive stars works for the first few Myr (Fig. [TJ af- 
ter star formation and might be attributed to observed anti- 
correlations of elements such as the Na-0 anti-correlation for 
star clusters with masses > 10 4 M@. 

Gas and metal return by supernovae occurs up to a 
few dozen Myr. As the ejecta by supernovae are much 
more energetic than winds by massive stars, star cluster 
masses must be high in order to keep a significant fraction 
of the material supplied by supernovae. Using en ergy argu- 
ments, iBaumgardt. Kroupa fc Parmentierl (|2008l ) derived a 
lower mass threshold of ~ 10 7 M© for star clusters to re- 
tain their residual gas despite multiple supernova events. 
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Figure 1. Sketch of different element and gas enrichment pro- 
cesses in massive star clusters. 



IWiinsch. Tenorio-Tagle. Palous~fc Silichl J 20081 ) show using 
2D-hydrodynamical simulations that also lower mass star 
clusters may be able to retain matter ejected by supernovae 
as a substantial fraction of such material can thermalise its 
high kinetic energy before escape from the star cluster. 

After all supernovae have exploded the massive AGB 
stars begin to evolve and are able to continuously replenish 
the gas reservoir in the star cluster and in its vicinity (e.g. 
iD'Antona fc Caloilliool ). 

After the epoch of massive energy feedback by the ra- 
diation and winds of massive ionising stars and supernovae, 
which prevents gas of the surrounding ISM from being ac- 
creted, the gas of the ISM is able to react to the gravitational 
potential of the new star cluster. It is expected that the more 
massive the star cluster is the stronger is the gravitational 
influence on the ISM by the new star cluster. In this paper 
we investigate this new possible scenario of gas accretion 
from the surrounding ISM by massive star clusters. 

Additionally, iFellhauer etHI (|2006h showed that if 
uj Cen formed in a dwarf galaxy then it could have cap- 
tured field stars during its formation from the underlying 
stellar field population contributing up to 40 per cent to the 
total final star cluster mass and showing the complex stellar 
population composition of the host galaxy. 

All existing explanations have in common that the ma- 
terial of further stellar generations has already assembled 
during cluster formation and thus their total mass does not 
increase with time. 

In contrast, Walc her et al. I (|2005h determined the 
masses of nuclear star clusters near the photometric cen- 
tre of bulge-less spiral galaxies to lie between 0.8 x 10 6 M© 
and 6 x 10 7 M© . As these nuclear clusters have luminosities 
by up to 2 orders of magnitude lar ger than the most lu - 
minous Milky Way globular clusters, IWalcher et all (|2005l ) 
suggested that these nuclear clusters grow by repeated ac- 
cretion of gas and show subsequent star formation. 

iBoker et al.l l|2004l) found that nuclear star clusters have 
dimensions (rm ~ 3.5 pc) comparable to massive globu- 
lar clusters, and some of these with masses >10 6 M© show 
evidence for multiple structural components and estimated 
periods between star formation episodes of about 0.5 Gyr 
l|Seth et al.ll2006l ). 
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3 MODEL 

From the hydrodynamical point of view the ISM as a self- 
gravitating gas is described by the three conservation equa- 
tions for mass, momentu m, and energy, the Poisson equation 
and the equation of state IjMartel. Evans II fc Shapiro 2006, 
e-g-)> 
dp 

at 



+ V • (pv) = 



<9v . VP „^ 
— + (v.V)v= V$, 

at p 

— + V • Ve = -PV • V H 

at P p 

V 2 $ = 4ttG(p - p) , 

p = fM , 



(i) 

(2) 

(3) 

(4) 
(5) 



where p is the mass density, P is the pressure, e is the specific 
energy, v is the velocity, $ is the gravitational potential, F 
is the radiative heading rate and A is the coolin g rate, p is 
the mean density. Following Mar tel et al.l (|2006l ) this term 
has to be included to prevent the overall collapse of the 
ISM and the term —p accounts for whichever process makes 
the ISM globally stable. If a compact massive star cluster 
is present within a homogeneous part of the ISM, then the 
potential in eq. 2 has to be split into a self-gravitating part 
and the external potential of the compact star cluster acting 
on the ISM. If the density fluctuations p—p are initially small 
then the total potential in the vicinity of the star cluster is 
mainly given by the star cluster itself. To further simplify 
the problem we assume that the star cluster has no relative 
velocity with respect to the embedding ISM, i.e. v = 0. With 
increasing relative velocity between the star cluster and the 
embedding ISM the effect of gas accretion is expected to 
become less dominant. Therefore, the stationary case gives 
a lower limit of a resulting mass threshold. 

As a final simplification only the static case can be con- 
sidered, i.e. all time derivatives vanish. This can be inter- 
preted such that the out-of-equilibrium ISM, for which the 
density distribution does not follow the static solution, at- 
tempts to reach this state. Thus, the static solution can be 
used to explore the gravitational influence of a compact mas- 
sive star cluster on the surrounding ISM and the full set of 
equations which are only numerically solvable reduces to the 
hydrostatic equation, 



-VP = -V*ci , 
P 

P = f(p, e) , 



(6) 
(7) 



where $ c i is the potential of the massive star cluster. 

The potential of the star cluster is assumed to be spher- 
ical. Then for a simple isothermal ideal gas, 



P = p^T. 
pm u 



(8) 



with a temperature T and a molecular mass pm u of the 
isothermal ISM, the static solution is given by 



p(r) = po e 



^ (*clM-*cl, ) 



(9) 



where po and & c \fi are the mass density and the star-cluster 
potential at a given reference radius, ro, introduced by the 



integration. The integration constants in eq. [9] are chosen 
such that the potential vanishes at infinity, i.e. for a large 
distance from the cluster centre, and the particle density is 
the particle density of the undisturbed ISM , 



p(oo) = pm u nisM , <&ci,o = $(oo) = , 



(10) 



where nisM is the particle density of the ISM. 

The potential created by the compact star cluster is 
descr ibed analytically by a Plummer potential (|Plummerl 

IS, 



$dW=-GM c ; 



+ &cl)~ 



(11) 



where M c \ is the mass of the star cluster and 6 c i is the 
Plummer parameter describing the compactness of the clus- 
ter potential and is e qual to the projected half-mass radius 
jrleggie fc Hudl2003h . 

The final expression of the static solution is 



n(r) = nisM e 



(12) 



The composition of the ISM is assumed to be primordial 
with X = 0.75 relative mass fraction of hydrogen and 
Y = 0.25 for helium. This corresponds to a mean molecular 
mass of p = A/ {AX + Y) = 1.23. In the following the solu- 
tions for ISM-star-cluster systems are plotted as a function 
of the star cluster mass and are parametrised by the ISM 
temperature and/or the Plummer parameter. If the ISM is 
enriched to for example Y = 0.3 then the mean molecular 
mass is p = 1.29, and the solutions correspond to the so- 
lution with the primordial composition but a star cluster 
mass reduced by 4.8 per cent. In the n-M c \ plot the curves 
then have to be shifted by 0.02 dex to the left. Thus the 
choice of a primordial ISM composition is appropriate for 
this analysis. 



4 THE 10 6 M Q -MASS THRESHOLD 

The Plummer parameter o f a star cluster is e qual to the 
projected half-mass radius (jHeggie fc Hutll2003l ). Assuming 
that the mass-to-light ratio of a massive star cluster is in- 
dependent of the radius then the projected half-mass radius 
is equal to the half-light radius. For Galactic globular clus- 
ters the median h alf-light radius is about 3.2 pc and ranges 
from 2 pc to 4 pc (Hasegan et al. 2005; Dabringhausen et al.l 
2008). The interstellar medium of the Galaxy consists of 
three main components l|McKeei Il995l ; iFerrierel 1200 ll ; ICoxl 
2005): a molecular and atomic cold component below 100 K 
but occupying only 1-2 per cent of the interstellar volume, 
a warm atomic and ionised component between 6 x 10 3 K 
and 10 4 K with a particle density of ~0. 1-1.0 cm -3 corre- 
sponding to the ionisation of hydrogen, and a hot component 
with temperatures with > 10 K and a particle density less 
than 10 -2 cm -3 fed by supernovae. The warm and the hot 
components have roughly the same volume filling factor of 
about 50 per cent. So it is most likely that star clusters are 
embedded in the warm or hot component of the ISM. 

The radial gas density profiles of three different star 
clusters {M c \ = 10 4 , 10 5 , 10 6 M@) with a Plummer parame- 
ter 6 c i=3.2 pc embedded in a warm ISM with a temperature 
of 8000 K and a particle density of 0.5 cm -3 are plotted in 
Fig. [2] It can be seen that the hydrostatic gas density within 
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Figure 2. Radial density profile of the static solution (eq. I Figure 4. Central density in dependence of the star cluster mass 

for three different star cluster masses (10 4 , 10 5 , 10 6 Mq) and for three different Plummer radii of 2, 3.2 and 4 pc embedded in 

a Plummer parameter b c \ = 3.2 pc. The embedding ISM has a an ISM with T = 8000 K and n = 0.5 cm -3 , 
temperature of 8000 K and a particle density of 0.5 cm -3 . 
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Figure 3. Central density in dependence of the star cluster mass 
with a Plummer radius of 3.2 pc for four different embedding 
interstellar media with a temperature of 6000 and 10 4 K and a 
density of 0.1 and 1 cm -3 . The horizontal line denotes the average 
density of matter in the solar interior. 



the Plummer radius does not vary much. We therefore use 
the central gas density, n c = n(r — 0), to characterise the 
static gas density in the inner part of a star cluster. 

In order to explore the influence of the cluster poten- 
tial on the embedding warm ISM we calculate the central 
density of the static solution (eq. I12|l in dependence of the 
star cluster mass with a Plummer radius of 3.2 pc shown 
in Fig [3] To cover the full range of the observed warm in- 
terstellar medium four different models are calculated with 
temperatures of 6000 K and 10 4 K and densities of 0.1 cm -3 
and 1 cm -3 corresponding to the boundary values of the ob- 
served warm ISM. For cluster masses smaller than 10 Mq 
the central particle density does not differ from the density 
of the undisturbed ISM. Between 10 5 and 10 6 Mq the cen- 
tral density starts to rise and becomes dramatically large 
above 10 6 Mq. 

The potential of massive star clusters with a mass larger 
than 10 6 Mq thus produces an instability in the warm ISM 



such that the ISM is expected to react with starting inflow 
towards the cluster centre. In other words: A self-gravitating 
ISM does not notice the existence of star clusters with a to- 
tal mass lower than 10 6 Mq , but star clusters more massive 
than 10 6 Mq become immediately attracting holes in the 
warm ISM. Note, how rapidly the instability rises with in- 
creasing cluster mass only depends on the temperature but 
not on the density of the ISM. Keeping the properties of the 
ISM constant (T = 8000 K and n = 0.5 cm" 3 ) the instabil- 
ity region also appears at cluster masses of about 10 6 Mq 
when varying the Plummer parameter between 2 and 4 pc 
(Fig.©. 

Present-day Milky Way globular clusters do not belong 
to the disk and are not embedded in a comoving warm ISM. 
But in early times when they formed they may have been 
embedded in a warm ISM, e.g. in a dwarf proto-galactic 
building block, ui Cen for example is believed to have been 
hosted in a dwarf galaxy disk before being accreted by the 
Milky Way. Because it was massive enough it could have 
accreted additional gas when the conditions were appropri- 
ate. Star formation events within the young ui Cen or in its 
vicinity are expected to have caused varying ISM properties. 
Thus it could have been placed alternatingly in the warm 
or hot ISM and the accretion history could have been fluc- 
tuating. Ionising OB stars in the young u) Cen would have 
prevented it from further accretion from the embedding ISM 
and star formation rested in uj Cen. Meanwhile star for- 
mation would have been ongoing in the ISM of the dwarf 
galaxy disk being continuously enriched with metals. After 
all ionising sources in u> Cen disappeared accretion would 
have restarted if the coomoving conditions were appropri- 
ate. When enough gas accumulated in ui Cen star formation 
would have restarted and the newly formed stars would have 
been more metal rich than the younger stellar populations. 
Thus, distinct populations with different metallicties are ex- 
pected in such a scenario in u Cen as observed. 

To underline the strength of this instability thresh- 
old for star cluster masses > 10 6 Mq the cooling time 
scale, t coo i, corresponding to the central gas density 
is calculated. We use the coolin g function provided 
in iKoeppenT Theis fc Hensierl (|l995h in the temperature 
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Plummer potential with mass M c i and Plummer parameter 
foci is 
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Figure 5. Cooling time-scale of the central static gas density in 
dependence on the star cluster mass for a Plummer parameter 
b c i = 3.2 pc. The embedding ISM has a temperature of 8000 K 
and a particle density of 0.5 cm~ 3 . 



regime from 100 to 10 4 K giving an estimate of the cool- 
ing time-scale of 



r C ooi/yr = 16400 (T/K) ' 5 (n/cm~ 3 ) 1 



(13) 



The resulting central cooling time-scale for the static solu- 
tion of an ISM with a temperature of 8000 K and a density 
of 0.5 cm -3 is plotted in dependence of the star cluster mass 
with a Plummer parameter of 3.2 pc (Fig. [5}. If the gas of 
the ISM starts to flow into the massive star cluster due to the 
instability caused by the cluster potential the increasing cen- 
tral gas density has a decreasing cooling time-scale support- 
ing the gas accretion. For comparison, a thermal instability 
leading to an e xpected coolin g flow of hot gas is known for 
galaxy clusters ([Fab ian 2003) and massive elliptical galaxies 
i|Kroupa fc Gilmordll994i r New star formation is inhibited 
by the ISM in massive elliptical galaxies being kept at a tem- 
perature of « 10 6 K due to the rando m stellar motions with 
veloc i ties of a few hundred km s -1 l|Mathews fc Brighentil 
2003:; IParriott fc Bregmanl liooih . The massive-star-cluster 
instability, however, leads to a cooling instability and the 
onset of star formation until the new OB stars reheat the 
cluster ISM. 

The argument can also be turned around. If warm ma- 
terial is released inside the cluster then the gas attempts 
to reach the hydrostatic solution. The warm gas of a star 
cluster less massive than 10 6 M@ will be distributed nearly 
uniformly in space, i.e. it escapes from the cluster, whereas 
a star cluster more massive than 10 6 Mq is able to keep its 
warm gas. 

It may be rashly argued that the reason for accretion 
or keeping of warm material is that the sound speed of the 
gas is smaller than the escape velocity from the star cluster. 
The sound velocity of an ideal gas is given by 



/ZTTlu S y (1 v K 



(14) 



For a monatomic gas, 7 = 5/3, the sound speed is about 
10.6 km/s (T = 10 4 K) or 8.2 km/s (T = 6000 K). The 
escape velocity at a distance of r = lb c \ from the centre of a 



'•- ■.-,/-^^= 92.7^(1 + 0"*, ( 15 ) 

VT+I foci S \l Mq fod 

The lowest threshold mass of 2.5 x 10 4 Mq is calculated in the 
case of the lowest temperature of the warm ISM (6000 K) 
and the central escape velocity (1 = 0). The largest threshold 
mass of 6.4 x 1O 4 M0 is calculated in the case of the highest 
temperature of the warm ISM (10 4 K) and an escape veloc- 
ity from the outer region of the cluster, e.g. the half mass 
radius (I = 1.305). Nevertheless, by comparing these two 
characteristic velocities describing the potential of the star 
cluster and the internal energy content of the gas thresh- 
old masses are obtained which are more than one order of 
magnitude smaller than the lower limit of masses of star 
clusters which multiple stellar population are observed in. 
In any case, the sound speed of a gas describes the velocity 
with which a pressure change propagates through space. No 
large scale motion of matter is involved in the propagation 
of sound, whereas the accretion of material is such a process. 

In any real situation of interest the gas will not be static 
and the interstellar medium will be turbulent and clumpy. 
Thus, it might be questionable how much is learned by com- 
puting static gas densities. Indeed, the calculated static so- 
lutions are not reached by the gas. The static solution de- 
scribes a stage of an equilibrium. In the case of the presence 
of a star cluster less massive than 1O 6 M0 this equilibrium is 
described by a nearly uniform density distribution. I.e. the 
warm ISM does not notice the existence of the star cluster. 
In the case of a star cluster more massive than 10 6 Mq the 
static solution is charactarised by a required central gas den- 
sity many orders of magnitude larger than the surrounding 
gas density. At the position of the massive star cluster a uni- 
formly distributed ISM is far from equilibrium. Thus, the gas 
tries to asymptotically reach this equilibrium by increasing 
its density at the position of the massive star cluster, i.e. the 
star cluster accretes. It will never reach the exact static gas 
density distribution as the increase of the central gas den- 
sity will result in a more efficient cooling and restarted star 
formation will stop further accretion by heating the infalling 
gas. 

Turbulences and inhomogeneities of the warm interstel- 
lar medium lead to condensation of material and the for- 
mation of cold molecular clouds at random locations in the 
galaxies. In this context massive compact star clusters with 
masses > IO^Mq stimulate this process at certain locations, 
i.e. the star clusters act as cloud condensation nuclei. 



5 THE AMOUNT OF ACCRETED MASS BY A 
MASSIVE STAR CLUSTER 

The threshold criterion derived here can be used to decide 
whether a star cluster is able to accrete gas from a warm 
medium or not. If a massive star cluster is able to accrete 
the question arises how much material can be accreted and 
if the amount of accreted material can account for the ob- 
served multiple stellar populations. Finding an accurate an- 
swer to this question requires extensive numerical simula- 
tion. However, a reasonable estimation can be done using 
the accretion rate expression for a star from its surrounding 
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medium (|Bondilll952h . The star has a potential of a point 
mass, whereas the potential of the star cluster is extended. 
But due to the spherical symmetry of the star cluster, the 
potential at a certain radius from the centre of the poten- 
tial is given by the enclosed mass within this radius placed 
at the origin of the potential (Newtons's second theorem). 
Therefore, outside the cluster radius the cluster potential is 
the same as if the cluster is assumed to be a point mass 
and the application of the Bondi-accretion is justified. The 
accretion rate, A, is given by 

A = 2ir (GM c i) 2 c sPoo , (16) 

where c s is the sound speed and poo is the mass density of the 
gas far from the cluster, i.e. the density of the surrounding 
ISM. With 

___ =2 . 9xl0 __ (17) 

and for a metallicity of \x = 1.23, a sound speed of c s = 
10.2 km/s and a particle density of nisM = 1 cm -3 the ac- 
cretion rate is about 3400 Mq /Myr. Thus on a time scale of 
a few 100 Myr corresponding to derived age gaps between 
the multiple stellar populations the possible amount of ac- 
creted material is of the order of 3.4 x 10 5 Mq for a 10 6 M© 
star cluster. The accreted material can account for the ob- 
served multiple stellar populations. 

We note in passing that the accretion may also be mod- 
ulated by the cluster orbit about its host galaxy. For example 
it may experience more significant accretion events when it 
passes through the ends of a galactic bar or through major 
spiral arms where the likelihood of it encountering cold gas 
is enhanced. 



6 ON THE COLD-GAS THRESHOLD FOR 
STAR CLUSTERS AND THE 
STAR-CLUSTER BIRTH INSTABILITY 

The criterion above which star cluster mass a density in- 
stability of the ISM can arise depends only on the ra- 
tio, n c /nisM, of the central and surrounding gas density in 
eq. 1121 The threshold expression is proportional to e MclllTbcl . 
Thus, it should be easier to accrete cold cloud material than 
warm material. But the volume filling factor of cold cloud 
material in galaxies is about 1-2 per cent and much smaller 
than the volume filling factor of the warm material (~ 50 
per cent). The possibility that a massive star cluster is long- 
term embedded in a cold molecular cloud is much smaller 
than the possibility that the star cluster is embedded in the 
warm ISM. If a star cluster is surrounded by cold dense 
material it is expected that accretion occurs for even less 
massive star clusters. As shown above the threshold is sur- 
passed for typical values of M c \ = 1O 6 M0, 6 c i = 3 pc , and 
T = 10 4 K. Star clusters should be able to accrete cold ma- 
terial with a temperature of about 100 K if the their mass- 
size ratio is larger than M c \/b c \ — 10 4 Mq/3 pc. For young 
embedded star clusters with b c \ = 0.3 pc (as the ONC for 
example) a total cluster mass of about lOOOM© is sufficient 
to accrete cold gas. 

It may be possible that this characterises the transi- 
tion from distributed star formation within giant molecular 
clouds in the form of loose groups hosting a few low-mass 



stars to star formation in well defined compact embedded 
star clusters: star formation in giant molecular clouds starts 
with distributed loose groups of a few low-mass stars. If 
regions in the low-mass-star forming giant molecular cloud 
surpasse the threshold expression, run-away accretion starts 
and a well defined star cluster appears. The ONC for exam- 
ple is a well defined young embedded star cluster within a 
molecular cloud showing many loose groups of few low-mass 
stars distributed over the molecular cloud. 



7 NUCLEAR STAR CLUSTERS 

Nuclear star clusters, which are embedded long-term in a 
warm comoving ISM and have masses less than 10 6 Mq , may 
not accrete additional material after their formation. But 
nuclear star clusters which are born or end up through merg- 
ers of smaller clusters more massive than 10 6 Mq should be 
able to accrete additional gas after their formation and show 
further star formatio n events. Therefore th ey are expected to 
grow as sug gested bv lWalcher et all (|2005l 'l. The star cluster 
mass region at 10 6 Mq should be underpopulated predomi- 
nantly for star clusters located in gas rich environments, i.e. 
the central regions of galaxies. 

Indeed, the ense mble of nuclear clusters observed by 
iGeorgiev et al.l (|2009h in gas rich dwarf irregulars have abso- 
lute V-band luminosities of —9.4 mag and less corresponding 
to stellar masses of w 1O 6 M and more. Compared with the 
population of blue and red globular clusters hosted in dlrrs, 
dEs and dSphs, which populate the distribution of V-band 
luminosities uniformly down to ~ —8.7, these nuclear star 
clusters seem to form a distinct population as they are sepa- 
rated from the rest of th e less massive GCs by a small gap in 
the luminosity sample (|Georgiev et alj|2009l . fig. 1). These 
nucear star clusters are candidates for repeated accretion 
events and so should contain complex stellar populations, 
as opposed to the clusters below the threshold which should 
not contain significant complex populations. 

It may be argued that nuclear star clusters are just fed 
with additional gas because they are located at the origin of 
the potential of the galaxy. But the location of nuclear star 
clusters and th e dynamical centers of disk galaxies do not al- 
ways coincide ( Mat thews fc Gallagher! |2002| ; IWalcher et al.l 
120051 ) . Recentlv lGeorgiev et all (|2009h found that seven of 
ten nuclear star clusters in dwarf irregular galaxies are off- 
centre by up to 480 pc. The overall process can be as follows. 

The central regions of galaxies are fed with gas due 
to the deep gravitational potential of the galaxy. For a nu- 
clear star cluster to be able to accrete gas from this continu- 
ously refilled central gas reservoir then depends only on the 
mass of the massive star cluster. Due to differential rotation 
the accreted material is expected to form a more rotation- 
ally than pressure supported new stellar system in the star 
cluster. Furthermore, younger stellar populations in these 
compact massive star clusters should be more elliptical and 
aligned with the galactic gas disk, whereas the initial pop- 
ulation should be more spheroidal if the formation of the 
star clusters occurred in a monolithic collapse. Indeed it has 
been found that older stellar populations in nuclear star clus- 
ters have a more spheroidal morphology and nuclear star 
clusters appear f lattened along the plane of the galaxy disk 
jSeth et alj|2008h . 
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Figure 6. Illustration of the half-light radius dynamical 
mass relation of pressur e supported stellar systems from 
iDabring hauscn et al. (2008). See text for details. 



iMilosavlievid l|2004) calculated the time-scale required 
for nuclear star clusters to have migrated from distant ec- 
centric orbits in the disk towards the central regions of their 
host galaxies. Such migration time-scales are too long and 
nuclear star clusters must have formed in situ in the cen- 
tres. However, if massive star clusters accrete gas from the 
long-term embedding ISM their masses grew. If the angular 
momentum of the star cluster is conserved their orbit must 
shrink and the inspiral may be accelerated. Places of the 
formation of such massive star cluster could have been for 
example the ma ssive clumps observed in chain galaxies at 
higher red shift (|Elmegreen fc Elmegreenll2006t ). 



8 OTHER MASS THRESHOLDS 

In the previous sections the influence of a few pc sized 
star cluster on the warm interstellar medium has been ex- 
plored. But such globular/nuclear star cluster type objects 
constitute only a subset of pressure supported stellar sys- 
tems, which do not have arbitrary dynamical properties 
but show distinct size-mass relations ( Hascgan et al. I [20051 ; 
iDabringhausen et al.ll2008l ; iForbes et al.ll2008f). 



Figure [6] is a sketch of figure 2 of Dabringhausen et al.l 

(2008) showing the half-light-radius-dynamical-mass rela- 
tions of pressure supported stellar systems. Here, the dy- 
namical mass, M, is the total mass within the optical extend 
of the object as obtained by solving the Jeans equations. 
Globular clusters do not show a mass-dependence of their 
half-light radius, whereas above about 10 6 Mq the half-light 
radius becomes mass dependent. The ultra compact dwarfs 
(UCDs), bulges of spiral galaxies and high- luminous ellipti- 
cals constitute a steep branch in the half-light -dynamical- 
mass plane. The radius-mass relation of these objects lying 
on the steep branch (SB) presented in IDabringhausen et al.l 
(2008) can be converted into a 6 c i — M relation, 



pc VlO 6 M e 



(18) 



The flat branch (FB) is build by dwarf spheroidals and 
low-luminous ellipticals. The b c \ — M of the flat branch is 

0.33 

(19) 
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Figure 7. Density ratio of the static g function of the total 

system mass of stellar systems lying on the fiat (FB) and steep 
(SB) branch (Fig.[HJ for the cases of the warm and hot interstellar 
medium. 



These two branches can be combined with the warm 
and the hot interstellar medium. The ratio of the central 
density and the density of the gas at infinity of the static 
solution feq, H2|l is plotted in Fig.[7]as a function of the total 
mass of the system. Systems with M > 10 6 Mq on the steep 
branch are able to keep their warm gas, whereas systems on 
the flat branch are able only to do so if their mass is 10 8 Mq 
or larger. This mass threshold coincides with the mass sep- 
arating the dwarf spheroidals and low-luminous ellipticals 
from each other. The hot interstellar medium can only be 
kept by systems more massive than about 10 Mq. This 
threshold mass corresponds to the transition region from 
low-luminous ellipticals on the flat branch to high-luminous 
spirals and bulges of disk galaxies on the steep branch. 

One may speculate if these coincidences of the mass 
thresholds obtained from the instability argument and the 
mass transitions observed in the mass-radius plane of pres- 
sure supported stellar systems has a direct physical origin. 



9 CONCLUSION 

We have shown that an ISM instability occurs for massive 
compact star clusters with masses of M c \ > 1O 6 M0. This 
instability may relate to episodic gas accretion from the em- 
bedding ISM and subsequent star formation. Such an ex- 
tended star formation history leads naturally to a spread in 
metallicity as observed in globular star clusters more mas- 
sive than « 10 6 Mq . 

This type of gas accretion can also account for the on- 
going star formation in massive nuclear star clusters which 
are still embedded in a comoving dense interstellar medium. 

The 'massive-star-cluster' instability adds to the well- 
known instabilities of the ISM, e.g. the Parker, Kelvin- 
Helmholz or Rayleigh- Taylor instability. 

Furthermore, for the combination of observed size-mass 
relations of pressure supported systems and different phases 
of the ISM the same analysis reveals other mass thresholds 
which coincide with transitions in the size-mass diagram 
of pressure supported stellar systems composed of globular 
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clusters, ultra compact dwarf galaxies, bulges of disk galax- 
ies, dwarf spheroidals and dwarf and large ellipticals. 

We thank Iskren Georgiev and J6rg Dabringhausen for 
helpful discussion on nuclear and globular star clusters. 
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